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APPENDIX I 

% Creating a network topology object 

network Jopo = topoCinit 1 ); % graphically place nodes on screen 

addlink(networkjopo); % graphically connect up nodes 

labelnames(networkjopo); % graphically label nodes 

save network_topo; % save network_topo for future use 



% Top level procedure to compute paths that optimize use of network capacity 
% inputs: 

% D = traffic demand matrix 

% (retrieved from predictions stored in TMS Statistics Repository) 
% network_topo = topology object defining the network topology 
% P = network policy information 

% (matrix of reserved capacity, which indicates links whose use 

% is administratively prohibited or which should not be 

% completely allocated) 

% outputs: 

% allocated_paths() = list of paths to set up, to TMS signalling system 



C = capacity(networkjopo); % retrieve network topology information 

C = C-P; 

saved_C = []; 
savedJSLA = []; 
assigned jpaths = []; 
round = 0; 

[SLA, S] = create_ordered_sla(D); 
F = SLA(1) 

for F - SLA', 

round = round +1; 
saved__C {round} = C; 
saved_SLA{round} = F; 

F % display the flow 

W = calc_weights( , calcweight2 , ,F,C); 
[dist, P] = floyd(W); 
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path = findpath(P,F.i ? F.j); 

assigned__paths {round} .path = path; 
5 assigned_paths {round} .flow = F; 

if (isempty(path)) 

fprintf(l, f no path for flow:\n'); F 

else 

10 C = compute_residual_capacity( f c - F.bw',path ? F,C); 

end 

end 

15 function [W] = calc_weights(func,F,C) 

% function [W] = calc_weights(func,F,C) 
% 

% Compute the weights by calling func on each elt of C 

% func must be of the form double ftmc(Flow F, Capacity_elt c, node i, node j) 
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func = fcnchk(func); 



fori = l:size(C,l) 

for j = l:size(C,2) 
25 W(i j) = feval(func ? F,C(i j),i j); 

end 

end 

function [w] = calcweight2(F,c ? ij) 
30 % function [w] = calcweight2(F ? c,i,j) 

% basic weight calc 

if (0 = c) 

w = inf; 
35 return; 
end 



% rule out paths that can't hack it 
40 if (F.bw>c) 

w = inf; 
return; 

end 
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w = 1 / (c - F.bw) ; % fill links with most capacity first 

function [C] = compute j:esidual_capacity(func ? path, F, C) 
% function [C] = computejresidual_capacity(func, path, F, C) 
% 

% Update capacity characteristics in C to reflect flow F being 

% allocated along path using function func 

% func should be of the form 

% C_element func(C_element c 5 Flow F) 

if(length(path) <= 1) 
return; 

end 

func = fcnchk(func 5 'c7F0; 

index = 1; 

src = path(index); 

index = index + 1 ; 

for index = index:length(path) 
dst = path(index); 

C(src,dst) = feval(func,C(src,dst),F); 
src = dst; 

end 

function [SLA, S] = create_ordered_sla(D) 

% function [SLA] = create_ordered_sla(D) 

% takes the demand matrix and returns a list of SLAs, 

% SLA of the form [ struct ; struct ; ... ] where struct is [BW ? 

% S of the form [ [BW, i, j] ; [BW, i, j] ; ■-.] 



S = []; 

for i = l:size(D ? l) 

for j = l:size(D 5 2) 

if(D(ij) — 0) 

S = [[D(ij)ij];S]; 

end 

end 

end 



36 



[Y, I] = sortrows(S,l); 

S = Y(size(Y,l):-l:l,:); % reverse order 

5 

SLA = struct( , bw^n^ml2cell(S(:J));i^num2cell(S(:,2));j',num2cell(S(:,3))); 
return; 

10 function [path] = findpath(P,i,j) 

% function [path] = fmdpath(P) 

% 
% 

15 path = []; 

if(i=j) 

path= [i]; 
return; 

20 end 

if(0 = P(ij)) 
path=[]; 

else 

25 path = [fmdpath(P,i,P(i,j)) j]; 

end 

function [D, P] = floyd(W) 
% function [D, P] = floyd(W) 
30 % given weights Wij, compute min dist Dij between node i to j 

% on shortest path from i to j J has immeadiate predecessor Pij 

n = size(W,l); 
if(n~=size(W,2)) 
35 error('Input W is not square??! !'); 

end 

D = W; 

40 P = repmat([l:n]',[l n]); 

P = P.*~isinf(W); 
P = P.*~eye(n); 

for k = 1 :n 
45 for i = 1 :n 
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l:n 

alt_path = D(i,k) + D(k,j); 
if(D(ij)>alt_path) 

D(i j) = alt_path; 

P(i j) = P(kj); 

end 
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APPENDIX II 

function addlink(TOPO) 
% addlink(TOPO) 
% 

% interactively add links to the TOPO 

update(TOPO); 
c_src= 1; 
c_dst = 2; 
c_bw = 3; 

figure(TOPO.curfig) 
while (1) 

fprintf(l,'\n\nHit Button 3 to end...\n\n'); 

% find coords and index i of src 
[xli yli button] = ginput(l); 
if (button ==3) break; end 

d = sqrt((TOP0.1ocs(:,l) - xli). A 2 + (TOP0.1ocs(:,2) - yli). A 2); 
[d,i] = min(d); 

xl = TOPO.locs(U); yl = TOP0.1ocs(i,2); 

% find coords and index j of dst 
[x2i y2i] = ginput(l); 

d = sqrt((TOP0.1ocs(:,l) - x2i). A 2 + (TOP0.1ocs(:,2) - y2i). A 2); 
[d,j] = min(d); 

x2 = TOP0.1ocs(j,l); y2 = TOP0.1ocs(j,2); 
hold on; 

lh = line([xl x2],[yl y2]/color7red'); 
cap = input('Enter capacity (in Mbps) > '); 

fprintf(l, About to create symetric %d Mbps link from node %d to node %d\n',cap,i,j); 

doit = input('Enter Y to confirm, N to reject, and B to change bandwidth (Y)> ','s'); 

if (isempty(doit)) doit = 'Y'; end 

if (doit == 'n' | doit == W) 
delete(lh); 
return; 
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end 

if (doit = V | doit == «') 

buf = sprintf('Enter capacity from %d to %d (in Mbps) > \i,j); 
cap_i_to J = input(buf); 

buf = sprintf('Enter capacity from %d to %d (in Mbps) > ' j 5 i); 
capj Jx>_i = input(buf); 

else 

ca PJLt°_J ~ ca p; 
capJ_to_i = cap; 

end 

%build the link records 
clear linkab linkba; 

linkab.src == i; 
linkab.dst=j; 
linkab.bw = cap_i_toJ; 
linkab.handle = lh; 

linkba. src = j; 
linkba. dst = i; 
linkba.bw = capJ_to_i; 
linkba.handle = lh; 



% now draw the actual link on the map 
delete(lh); 

lh = drawlink(TOPO ? linkab); 

% now store the link info 

TOPO.links = [TOPO.links ; linkab ; linkba ]; 

TOPO.linkarray = [TOPO.linkarray ; [ i j capjjoj] ; [ j i capJ_to_i 

end % of while loop 
assignin('caller r ? inputname(l) ; TOPO); 



function [C, portmap] = capacity(TOPO) 
% [C, portmap] = capacity(TOPO) 
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% portmap maps indices of C to elts of nodes(TOPO) 
% [node dir] where 

% node is index of elt in nodes(TOPO) 

% dir is 1 if data enters here, -1 if data leaves here 

5 

numnodes = length(TOPO. links) * 2; 

C = zeros(numnodes,numnodes); 

10 curnode ^O; 

portmap = []; 

for i = 1 :length(TOP0.1inks) 
link-TOPO.links(i); 
curnode = curnode + 1 ; 
15 portmap(curnode ? :) = [link.src -1]; 

curnode = curnode + 1 ; 
portmap(curnode,:) = [link.dst 1]; 

C(curnode- 1 ,curnode) = link.bw; 

^tO end 

S P c_node == 1; 

,h c_dir = 2; 

ii 25 for i = 1 :length(TOPO.nodes) 

O ins = find(portmap(: ? c_node) == i & portmap(: 5 c_dir) == 1); 

M; outs = find(portmap(:,c_node) == i & portmap(:,c_dir) == -1); 

!tl forj^ins 

f^O for k = outs 

CG,k) = inf; 

end 

end 

end 

35 function [a, b ? c] = debug(t) 

update(t); 
fieldnames(t) 

40 

a = t.nodes 
b = tlocs 
c = t. links 

function display(TOPO) 
45 % DISPLAY a topo object 
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% a link is a unidirectional, so the value is probably twice what you want 

fprintf('[TOPO object: %d nodes %d links]\n\... 

length(TOPO.nodes) ? length(TOP0.1inks)); 
function draw(TOPO) 
% draw(topo) 
% 

% draw the topology figure in a new window 

TOPO.curJIg = figure; 

axis(TOPO.axis); 

axis equal; 

axis manual; 

box on; 

hold on; 

for i = 1 :length(TOPO.nodes) 

run = plot(TOPO.nodes{i}.loc(l) ? TOPO.nodes{i}.loc(2);ob , ); 
TOPO.nodes{i}.mark_handle = nm; 
if(isfield(TOPO.nodes{i} 5 , nameloc')) 

TOPO.nodes{i}.nameloc(3) = text(TOPO.nodes{i}.nameloc(l),. 

TOPO . nodes { i } ,nameloc(2),TOPO .nodes { i } .name 

end 

end 

% yes, this draws the same link twice, fix it if it matters -dam 11/21 

TOP0.1inkarray = []; 

for i = 1 :length(TOP0.1inks) 

TOP0.1inks(i),handle = drawlink(TOP0 5 TOPO,links(i)); 

TOPO.linkarray = [TOPOlinkarray ; ... 

[ TOPO,links(i).src TOP0.1inks(i).dst TOP0.1inks(i).bw]]; 

end 

assignin('caller T 5 inputname(l) ? TOPO); 
function ex(t) 

t.nodes 

function labelnames(TOPO) 
% function labelnames(TOPO) 
% make it easy to label the nodes 

for i = 1 :length(TOPO.nodes) 

fprintf(Tlace label for node %d ?, %s 1, \n , J i J char(TOPO.no<ies{i}.name)); 
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origcolor = get(TOPO.nodes{i}.mark_handle 5 ! color r ); 
set(TOPO.nodes{i}.mark_handle;color' 5 [l 0 0]); 

if (isfield(TOPO.nodes{i};nameloc')) 

good_x = TOPO.nodes{i}.nameloc(l); 
goodly = TOPO.nodes{i}.nameloc(2); 

end 
th = []; 
while (1) 

fprintf('Button 1 to (re)place text, Button 3 to accept\n'); 

[x,y,button] = ginput(l); 

if (3 = button) break; end 

if (~isempty(th)) delete(th); end 

th - text(x,y,TOPO.nodes{i}.name); 

good_x = x; goodly = y; 

end 

TOPO.nodes{i}.nameloc = [good_x, good_y, th]; 
set(TOPO.nodes{i}.mark_handle,'color',origcolor); 

end 

assignin( t caller l ; inputname(l) 5 TOPO);function names(TOPO) 
% NAMES the list of names of the nodes in the topo 

fprintf('Node\t\tName\n'); 
for i = 1 :size(TOPO.names,l) 

fprintfC%d\t\t%s\n' ? i 5 TOPO.names{i}); 

end 

function [node] = nodes(TOPO) 

% function [node] = nodes(TOPO) 

% returns a cell array describing nodes in the TOPO 

node = TOPO.nodes; 

function [TOPO] - topo(TOPO) 

%[TOPO]-topo(TOPO) 

%% if input TOPO is 'init 1 , create a new topology 

% 

% newtopo = topoCinit'); 

% 

% else add new nodes to TOPO 

% 

% nodes is a array of structs, one per node 

% link is a array of structs, one per link 

% a link is a unidirectional item, so there are probably twice 

% as many links as you'd expect. 
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if ( Mr 8^^ opo(TOp0) or t opo("init") - not enough args') 
end 

if (ischar(TOPO) & TOPO == 'initf) 
clear TOPO 

TOPO.nodes = []; 
TOPO. links = []; 

TOPO.capacity = []; % now computed as needed 
TOPO.locs = [] ; % internal cache 
TOPO.linkarray - []; % internal cache 

f = figure; 
axis([0 75 0 50 J); 
TOPO.axis = axis; 
TOPO.cur_fig = f; 
axis equal 
axis manual 
box on 
hold 



else 
end 



figure(TOPO.cur_fig); 



nodecount = length(TOPO. nodes); 

while (1) 

clear nodeinfo; 

fprintf(l,'\n\nHit Button 3 to stop\n\n'); 

[x y but] = ginput(l); 

if (but = 3) break; end 

x = floor(x); y = floor(y); 

nm = plot(x,y,'ob'); 

name = input('Enter name > ','s'); 

nodeinfo. loc = [ x y]; 
nodeinfo. mark_handle = nm; 
nodeinfo.name = cellstr(name); 
nodecount = nodecount + 1 ; 
TOPO.nodes{nodecount} = nodeinfo; 

end 
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if (topo' ~= class(TOPO)) 

TOPO = class(TOPO,'topo'); 

end 

5 

if (nargout = 0) 

assignin('caller' ,inputname( 1 ),TOPO) ; 

end 

function lh = drawlink(TOPO, link) 
10 % assumes TOPO.linkarray is already valid, and draws the position of 

% link line based on the number of links already present in linkarray 

c_src = 1 
c_dst = 2 
15 c_bw = 3 

i = link.src; 
% j = link.dst; 

Jo xl =TOPO.nodes{i}.loc(l) 

2 yl =TOPO.nodes{i}.loc(2) 
C x2 = TOPO.nodes{j}.loc(l) 
§ y2 = TOPO.nodes{j}.loc(2) 

L25 

=j if (isempty(TOPO.linkarray)) 
Z numjinks = 0; 

|f else 

3 numjinks = sum(TOP0.1inkarray(:,c_src) = i & TOP0.1inkarray(:,c_dst) == j); 
130 end 

pattern=[0 1 -12-2 3-3] * .3; 

if(abs(xl -x2)>abs(yl -y2)) 
35 delta_x = 0; 

delta_y = pattern(num_links +1); 

else 

delta_x = pattern(num_links +1); 
delta_y = 0; 

40 end 

lh = line([xl x2] + delta_x, [yl y2] + delta_y, 'color', 'black'); 
function update(TOPO) 
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clear TOPO.locs; 

for i = l:length(TOPO.nodes) 

TOP0.1ocs(i,:) - TOPO.nodes{i}.loc 

end 

clear TOPO.linkarray; 

for i = l:length(TOP0.1inks) 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ TOPO.linksQ.src TOP0.1inks(i).dst TOP0.1inks(i).bw]]; 

end 

% these are here to be cut and pasted into other functions as needed 

% there doesn't seem to be a good way to pass them around in another fashion 

% (using assigning('caller\..) to force their definition sounds like asking 

% for trouble 'cause you'll overwrite another definition of them...) 

c_src - 1 ; 

c_dst = 2; 

c_bw = 3; 

assignin( f caller , ,inputname(l),TOPO); 



